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1/^ ' Abstract 

(N ■ 



Motivated by certain mathematical models for Micro-Electro- Mechanical Systems (MEMS), we give upper 
and lower L°° estimates for the minimal solutions of nonlinear eigenvalue problems of the form —Aw — 
\f{x)F{u) on a smooth bounded domain in R^. We are mainly interested in the pull-in distance, that is 
1^ ' the norm of the extremal solution u* and how it depends on the geometry of the domain, the dimension 

I of the space, and the so-called permittivity profile f. In particular, our results provide mathematical proofs 

for various observed phenomena, as well as rigorous derivations for several estimates obtained numerically 
by Pelesko [T7], Guo-Pan-Ward fTJ" and others in the case of the MEMS non-linearity F{u) = (^i\y2 and 
for power-law permittivity profiles f{x) — |a;|". 

> \ 

^ ■ 1 Introduction 

' We examine problems of the form 

cn 
o . 

^ . 1 u = ondn, 

O 

where £7 is a bounded domain in M^, < A, / is a nonnegative nonzero bounded Holder continuous function, 
usually dubbed as the permittivity profile, and F is a smooth, increasing, convex nonlinearity on its domain 
€ Dp C M, such that F{0) = 1 and which blows up at the endpoint of its domain. We shall concentrate 
' on the two cases where either F is superlinear and its domain is Dp := [0, -l-oo) in which case F is said to 

be a regular nonlinearity, or when Dp := [0, 1) and liuiuy'i F{u) = +oo in which case, we say that _F is a 
singular non-linearity. Typical regular nonlinearities are F{u) = e" or F{u) = (1 + for p > 1, while 
singular nonlinearities include F[u) = {1 — u)~^ for p > 0. 

We say that a solution u of (-Pa,/) is classical provided 111*11^°° < oo (resp., < 1) if is a regular 

(resp., singular) nonlinearity. Note that by elliptic regularity theory, this is equivalent to saying that a 
classical solution is in C^'" for some a > 0. 

We shall also need to consider HQ — weak solutions of (Pa,/) which arc those u in iJo(ri) such that 

/j^ VuV(t)dx = /j-^ \fF{u)(j)dx for aU e Hl{n). (1) 



-Au = Xf(x)F(u) in n. 
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It is by now well-known that - regardless whether is a regular or singular nonlinearity - there exists 
an extremal parameter A* € (0, -|-cxd) depending on Q, f and N, and which can be defined as 



A*(r2, /) := sup{A > : (-Pa,/) has a classical solution}, 

such that (-Pa,/) has a minimal classical solution ua for every A G (0, A*), and no weak solution for A > A*. 
By a "minimal solution" u, we mean one such that any other solution v of (Pa,/) satisfies v > u a.e. in Q. 
One can then also show that A i-^ u\{x) is increasing on (0, A*) for each x & ft. This allows us to define the 
extremal solution by 

u*(x) :— Van u\(x), 
xyx* 

which then can be shown to be the unique (weak) solution of (Pa*,/). 

We shall also need the notion of stability. Given a weak solution u of (Pa,/), we say that u is stable (resp., 
semi-stable) provided /ii(A,M) > 0, (resp., /zi(A, u) > 0) where 

fii{X,u) := inf ( / dWI' - Xf{x)F'{uyj^)dx : ^ € H^{n), / = 
iJn Jn 

Under our assumptions on the nonlinearity F, and whether it is regular or singular, one can show that for 
all < A < A* the minimal solution u\ is stable, and consequently that u* is semi-stable. If in addition, u* 
is a classical solution of (P\*,/), then necessarily /ii(A*,u*) = 0, since otherwise one could use the Implicit 
Function Theorem, in a suitable function space, to obtain solutions to (-Pa,/) for A > A*, which would be a 
contradiction. On the other hand, one has the following useful result, which was proved by Brezis-Vasquez 
[1] for regular nonlinearities, and by Ghoussoub-Guo JlOj in the case of singular nonlinearities with general 
permittivity profiles. 

Proposition 1.1. A semi-stable H^{Vt)~weak solution of {P\j) that is not a classical solution can only 
occur at A* , in which case it must be equal to the extremal solution u* . 

The question of the regularity of the extremal solution has attracted a lot of attention in the last decade. 
For general regular nonlinearities the extremal solution is classical provided one of the following holds: 

• 57 is contained in with < 3 (Nedev, see |16)). 

• 17 is a ball in with iV < 9 (Cabre and Capella, see P]). 

The second result is optimal after one considers F{u) — e" on the unit ball in M^". It is an open question as to 
whether for 4 < ./V < 9, there is a regular nonlinearity F and a domain 17 C on which the corresponding 
extremal solution is unbounded. In the case of the MEMS model, where F{u) = (1 — it)~^, it is known that 
the extremal solution is classical provided N < 7 and that this result is optimal (see [10]). On the other 
hand, for any dimension N > 2, there exists a singular nonlinearity, namely F{u) = (1 — u)^^ for some 
p :— p{N) > 0, such that the corresponding extremal is not classical (see Chapter 3 of Ff]). 

In this paper, we are mostly interested in the quantitative aspects of the regularity of the extremal 
solution u* , which were initially motivated by the equation 



-Au ^ in f7, 

^ (Ma,/) 



on dn. 



In dimension N = 2 this equation models a simple Micro-Electromechanical-Systems MEMS device, which 
roughly consists of a dielectric elastic membrane that is attached to the boundary of 17, and whose upper 
surface has a thin conducting film. At a distance of 1 above the undeflected membrane sits a grounded plate, 
i.e., a plate held at zero voltage. When a voltage ^ > is applied to the thin film of the membrane, it deflects 
towards the ground plate. After various physical limits of the parameters involved, a dimensional argument 
and a simplification, ones arrives at (Ma,/) for the steady state of the membrane. Here A is proportional 



2 



to the applied voltage V and the permittivity profile f{x) allows for varying dielectric properties of the 
membrane. 

As seen above, one expects the extremal solution u* in small dimension N to be bounded away from 1, 
hence to be a classical solution. Since the parameter A* corresponds to the critical voltage beyond which 
there is a snap-through, and since u* is the optimal deflection of the membrane, it is therefore important 
for the design of MEMS devices to know how the critical voltage A* and the pull-in- distance - defined as 
||m*||2,oo - depend on the geometry of the membrane and on the permittivity profile. Several analytical and 
numerical estimates on A* have been derived by Pelesko [T7], Guo-Pan-Ward [13], Guo-Ghoussoub [TO] and 
others in the case of the MEMS non-linearity F{u) = (^i\yi ■ On the other hand, only numerical estimates 
have been obtained for the pull-in distance in the case of power-law (resp., exponential) permittivity profiles 
f{x) = (resp., f{x) = e"^). In this paper, we shall see that one can give rigorous proofs and estimates 
for phenomena, which so far have only been observed numerically by various authors. We shall also include 
corresponding results for general ~ not necessarily MEMS-type - nonlinearities. 

Here is a brief description of the paper. In section 2, we give upper estimates on the pull-in voltage A* (O, /) 
in fairly general situations, which will in turn yield lower bounds on ||M*||ioo. What is remarkable here is 
that the estimates - which are valid for general nonlinearities - turn out to only depend on the permittivity 
profiles and not on the domain, nor on the dimension. Actually, they also apply to any reasonable uniformly 
elliptic operator. 

In section 3, we give upper estimates on ||u*||l=° which are computationally friendly. Just as in the proof 
of the regularity of u* in low dimensions, we use the energy estimates on the minimal solutions coupled with 
LP to L°" Sobolev-type constants related to corresponding linear equations. While the result is satisfactory 
for exponential nonlinearity, it is not so for the MEMS model, which led us to reconsider this nonlinearity 
in the case of the ball where more precise to Holder estimates can be used. We stress here that we are 
not interested in optimal upper estimates but rather estimates which, if given a specific domain and a 
nonlinearity _F, one can easily obtain some numerical parameters by plotting a function of a single variable 
- possibly - using a Computer Algebra System. 

Section 4 was motivated by an intriguing phenomena observed numerically by Guo-Pan-Ward [13| , namely 
that on a two dimensional disc, the pull-in distance does not depend on the power of the permittivity profile 
f{x) = |a;|". We prove that this is indeed the case by a simple scaling argument which relates the problem 
(-fA,|x|o) on the unit ball of to (^^,1) (which for simplicity we denote by (Px)) on a ball in a fractional 
dimension N{a). (Note that when / is radial and Vl is the unit ball in R^, all stable solutions of {P\.f) are 
then radial and hence we can examine the problem in fractional dimensions). One can then easily transfer 
many results established for (Px) to (P>^ |^|q). This observation, combined for example with the results of 
Cabre and Cappella [2], leads to new regularity results for the extremal solution associated with (-Pa.|x|°)- 

In section 5, we study the asymptotics in A, and we obtain upper and lower pointwise bounds on the 
minimal solutions ma, in the case where u* is singular. The upper estimates are valid on arbitrary domains 
and we restrict ourselves to radial domains for the lower estimates since more explicit bounds can then be 
found. For that, we show that A 1— > is actually convex, and we exploit the fact that both u* and ^ua|a=a* 
are explicitly known in the case where 17 is a ball and u* is singular. 

We now list our main notation. For a nonlinearity F, we denote by a^? the upper bound of the domain 
Dp, which means that 00 if is regular, and ap ■= 1 if is singular, in such a way that Dp := [0, of). 

We shall also associate to F the numbers 

Bp := sup ^ and Cp -.^ C (2) 

Te{0,ap) 

The ball of radius R centred at xq in R^ will be denoted by Bii{xo). li xq = then we omit a;o and if i? = 1 
then we just write B. Given a set in R^ we let |0| denote its A^-dimensional Lebesgue measure, while o^jv 
denotes the volume of the unit ball B in R^. The conjugate index of p will be denoted by p' in such a way 
that - + ^ = 1. For a radial function u we write u{r) = u(|a;|). The first eigenvalue of —A in HQ{n) will 
be denoted by Ai(f2) and the corresponding positive eigenfunction will be (f)n, assuming the normalization 
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2 Lower estimates for the L^— norm of the extremal solution 

This section is devoted to the proof of the foUowing result. 

Theorem 2.1. Suppose F is either a regular or singular nonlinearity and that u* is the extremal solution 
of {P\j), which we assume to be classical. Then, 

where we define {F')^^{z) — for z < F'{0). 

Before proceeding with the proof, we give some apphcations. 

Corollary 2.1. Suppose f is a non-negative bounded Holder continuous permittivity profile and that the 
extremal solution u* of (Px.f) on a bounded domain Q is regular. 

1. IfF{u) = (jtTijF; P > 0, then 

[p+lVnifo// \p+l j^fcjyndxl J 

In particular, when the permittivity / = 1, then for any dimension 1 < N < 2 + + 4y/^^, and 
any bounded domain Q, C M^, we have 

\\u*\\l^ > (5) 
p+l 

2. IfF{u) = (m + l)P, p>l, then 

p ( info fP~ ^ In f4>ndx \ ^ 



u > max 



p - 1 Vsupo/ 



) \ p supo / / J 



In particular, when / = 1, then for any dimension 1 < N < 10, and any bounded domain fl C M^, we 
have > 

3. IfF{u) = e", then 



>max|l + log( ^"^"{. ), logf 



In particular, when the permittivity f = 1, then for any dimension 1 < N < 9, and any bounded 
domain fi C M.^ , we have \\u*\\Lca > 1. 

Remark 2.1. Note that the lower bounds (when / = 1) are independent of the domain. It is also fairly 
easy to adapt the proof below to show that they are not particularly exclusive to the Laplacian — A. Indeed, 
the same lower bounds can be obtained if we replace it by any operator of the form L(u) := — div(A(x)VM) 
where A{x) is a symmetric uniformly positive definite N x N matrix defined in Q. 
Moreover, the same arguments show that the extremal solution associated with 

A^u = \F{u) on n, (8) 

also satisfies the same lower bound, where for general domains ft we restrict our attention to the Navier 
boundary conditions: u ~ Au = on dfl, while in the case of Q being a ball we can use the Dirichlct boundary 
conditions: u = d,^u = on dB. For recent advances on fourth order nonlinear eigenvalue problems, we refer 
to [3], [5, and [6j. 
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The proof of Theorem l2.1l foUows immediately from the combination of the following two propositions. 
The first provides upper estimates on A*, in terms of F, f2 and /. 

Proposition 2.1. Suppose F is either a regular or singular nonlinearity. Then 

A-(n,/)<A,(n)™{^,^}, ,8) 

where Bp and Cp are given in (0). 

Proof. Supposing it is a classical solution of {P\j), we multiply both sides of the equation by 0o and integrate 
to obtain 



{XF{u)f - Xi{n)u) (j)ndx = 0. 



Since 0n > we must have 



A < Ai(f})sup-— — - < — — sup -— = 

n jF{u) mt(i J zgDf F[z) mffi / 

For the second bound, multiply (Pa./) by and integrate to obtain 

l-Au) — 
n ^ ( 

2 



Xf(t)ndx = / {-l\u)——^dx 
r [uj 



F{u) F{u 
Ai(r!) I (j>n\ I ITT-^dA 



n 



Fir) 



< Ai(17)Cf, 

after recalling the normalization of (pn- □ 
Proposition 2.2. Suppose u* is the extremal solution of (Pa,/) which we assume to be classical. Then 

Xi{n)< \*\\fF'{u*)\\L^. (10) 

Proof. If w is a classical solution of (Pa,/) with Ai(il) > X\\ f F' (u)\\ l°° , then the variational formulation of 
the first eigenvalue Ai(f7) of the Laplacian yields 

\W(j>\^dx> Xi{n) / (p^dx > X\\fF'{u)\\L'^ / (f>^dx>X / fF'{u)(t)^dx, 
Jn Jn Jn 

which means that u is then a stable solution of (Pa,/)- 

Assuming now that u* is a classical solution then, as mentioned in the introduction, we necessarily have 
that fj,i(u*) = 0. Using the bifurcation theorem of Crandall-Rabinowitz one can then obtain a second 
branch of solutions U\ to (Pa,/) for A in a small interval (A* —e, A*). Moreover these solutions are unstable in 
the sense that ^i(A, U\) < 0. It then follows from the above that Xi{n) < A||/P'(?7a)||l°°- Sending A Z' A* 
gives the desired result. □ 
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3 Upper estimates for the L^— norm of the extremal solution 

In this section we look for upper estimates on the extremal solution u* associated with (Pa), where F is one 
of the three linearities considered in Corollary 12. 11 and where we take f{x) — 1 for simplicity. The methods 
consist of combining the energy estimates - which are critical in showing that the extremal solution is regular 
in low dimension - with various L°° and Holder estimates for linear equations. 

The following simple observation can be useful when looking for upper estimates. 

Observation 3.1. Suppose u* is the extremal solution associated with {P\) in il with extremal parameter 
A*. Then the extremal solution associated with {P\) in the domain Qp :— pCl (where p > 0) is given by 
v*{x) := u*{^) with extremal parameter A* (lip) = 

3.1 Upper estimates on general domains 

We begin with the case of exponential nonlinearities. 

Theorem 3.1. Suppose F{u) = e", is a bounded domain in and u* is the extremal solution associated 
with (Pa). 



1. //3 < TV < 9, then 

\\u\\l'^<—^ — ^( — J ' (11) 

where 

2t . 2t / 4 \i iV-2 

— — <t <2 



e{N -2)\ujnJ 

:= inf \ ( )^ f ^) ' : 

^ 1 \4:t + 2-N' \2-tJ 

2. IfncBiC M2, then 

where we define for p > I and < R < 1, 

A{p,R) := / {~\og{r))Prdr. 

Using a computer algebra system one can evaluate the constants and obtain: 

Ps = 1.9915, Pi 2.2324, /Jg 2.6689, Pe = 3.42269, 

/37 = 4.81191, /Jg = 7.9408166, ^9 = 19.0031. 

Note that one can combine this upper estimate with the previous lower estimate on u* to obtain the following 
lower bound on the first eigenvalue of the Laplacian on a bounded domain in M.^ whenever 3 < TV < 9: 

We also consider the case of a MEMS nonlinearity. 

Theorem 3.2. Suppose F{u) = (1 — u)^^, fl is a bounded domain in M.^ and u* is the extremal solution 
associated with {M\) infl. If3<N<7, then 



< 1 
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where 



, \ IQN^ f 2t \^/4(2t+l)\' 3(A/'-2) „ ^ 

Remeirk 3.1. f/sm^ a similar approach one can show that if u* is the extremal solution associated with 
(Pa); in the case where F{u) = {u + ly, p> 1, and N = 3 or N = A then 

^ (p-lf-Ui(0)/3iv,p /|n|^^ 
II" 11^-^ pP{N-2) ' 

where 



I3n,p = mf <^ ^ ^ , ^ ^ ^ : max{t;, ijv,p} < t < i + 



anrf where 

r~n + i~o P 1 

tp ■■=P- VP -P, tl :=p+ - P, tN,p ~~^~2'^2' 

We have omitted N = 2 just for simplicity. To obtain estimates for N < 10 one has to perform a bootstrap 
argument or restrict the range of values for p. 

For proving the above theorems we shall need the following easy lemmas. 

Lemma 3.1. Let Cl be a smooth bounded domain in M^. 

1- If N >3 and t> ^, then for all x € il, 



jdy ) < 



L\y-x\'^^-^)^' V " (2r-iV) 
2. If N = 2 andfl C Bi c K^, then for all xeQ, 



log(|y - x\)r'dy^ ^ < (2.)^a(-^, i^) . 



We now obtain L°° bounds on linear equations. 

Lemma 3.2. Suppose —Aw = g{x) > m with u = on dil where Q a bounded domain in and g is 
smooth. 

1. IfN> 3, then for all r > f , 

yk^(r-l)^|0| 



1 , „ , 2 1 

Noll rrlr — 1 ) T 

■U TOO < 



Nt{N-2)uj^ (2r-iV)- 

2. If N = 2 andn<zBx, then for all r > 1, 
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Proof. In both cases, we let v{x) denote the Newtonian potential of g, i.e., 

^("^^ ■= N{N -2)u;m L 
for AT > 3 and 

v{x) '■= ^ J^i- logdy - x\))g{y)dy, 

for N = 2. Since < u{x) < v{x) in f2, it suffices to show the desired L°° estimate on v. To do this, one 
uses (for A'^ > 3) Holder's inequality to write 

and then use the integral estimate in the previous lemma. □ 

We now derive the energy estimates for stable solutions. 
Lemma 3.3. Suppose u is a classical semi-stable solution of {P\). 
1. If F{u) = e", then for all < t < 2, we have 

4 



2. If F{u) = (1 - then for all0<t<2 + we have 



in 

- ^2tu 



Proof. 1) Using the test function ^ := e*" — 1, where < t < 2, in the stability conditions gives 

Jn 

Now testing (P\) on </> = e^*" — 1 and rearranging, gives 
Comparing the last two inequalities and dropping some positive terms gives 

and after an application of Holder's inequality on the right one obtains 

||e"|k.*+i <--^|J7|^. (15) 
(2 — t)t 

2) Take ?/> := (1 — w)"^ — 1, (6 := (1 — ?i)~2*-i _ i ^nd proceed as in 1) by putting tp into the stability 
condition and testing (Ma) on (p. We obtain 



(i2 2t+l) L{l-ur^+^- t^ Ja 



after dropping a couple of positive terms. Holder's inequality then yields 



/2_ 1 ^ 


1 


t 


U2 2t + lJ 


1-u 


Z,2t+3 



<l|n|^. (16) 



□ 
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We now combine the energy estimates with the hnear estimates to obtain upper estimates on u* . 



2t + 1 along with the estimate 



Proof of Theorem I3.lt Use Lemma [3.21 with g{x) :— A*e" and r 
A* < to arrive at an estimate of the form 

— e 

\\u*\\l^ < C(t,iV,|r!|)^i^||e"*||i2t+i, 
e 

where C{t, N,\fl\) is provided by Lemma [3.21 Now replace the L^*— norm on the right using the energy 
estimates from Lemma [5751 to arrive at the desired result. The restrictions on t are a result of the restrictions 
on T in the linear estimates along with the restrictions on t from the energy estimates. □ 

Proof of Theorem 13. 2t Let il denote a bounded domain in M.^ where 3 < < 7 and let u* denote the 
extremal solution associated with (Mx) in O. Since the reasoning works for any log-convex nonlinearity F 
(i.e., u 1-^ log{F{u)) is convex), we define v := log(F(u*)), and so 

-Av = - — 
dv? 

with w = on d^l. Since F is log convex, the first term on the right is negative. We now define w by 

-Aw = \*F'{u) in 
w = Q on 



2^og{F{u))\ 



\*F'(u 



in ri, 



and so < v{x) < w{x) a.e. in f2 by the maximum principle. Using the linear estimates from Lemma 
with g{x) := X*F'{u*) one has 



log 



1 



from Lemma 



which will give the desired conclusion. 



1 



1 - u* 



^2,U- = IklU- < M\l^ < Cr\*\\F\u*)\\L^ =CrX* 

Taking now r = ^ + l>-^,we can then replace the L'^— norm on the right by using the energy estimates 



□ 



3.2 Upper estimates on radial domains 

While the upper estimate on general domains obtained in the last subsection is quite satisfactory for the 
exponential nonlinearity, it is not so for the case of the MEMS nonlinearity. Indeed, using Maple one sees 
that if :— (0, 1)^ the unit cube in (and so Xi{fl) — 371^), Formula (fTij) would then give that 

\\u*\\l^ < .993..., (17) 

which is clearly not a very good upper estimate. This is mainly due to the fact that we drop a potentially 
large term in the proof of Theorem 3.2, when we replaced w by it; in order to apply the linear estimate of 
Lemma 3.2. Note that this was not needed for the exponential nonlinearity in the proof of Theorem 3.1. 

In this section we examine radial domains, where better results are available on u* , at least in the case of 
F{u) = (1 — w)~2. One can also examine the exponential nonlinearity using this approach but we won't do 
this since the last section seems to give satisfactory results. For simplicity, we shall also restrict our attention 
to the case of / = 1. The main difference is that we use here Holder estimates on linear equations versus 
the L°° estimates of the last subsection. 

For the remainder of this section we assume that H. is the unit ball B in and F{u) ~ (1 — u)^^. We 
define the following parameter: 



Lemma 3.4. Let u denote a smooth radially decreasing solution of —Au = g{r) > Q in the unit ball B of 
M^. //max{l, -y} < r < oo, then one has the estimate: 

u[Q) > u{R) > u{0) - 7(r.iV)|g|U. j^2-f for all Re (0,1). (18) 

Proof. When = 1, we integrate the equation between and r, and apply Holder's inequality to obtain 

1 

—u'{r) < ^sILlU—, Now integrate both terms between and R, and use again Holder's inequality to obtain 
the desired result. 

When N > 2, we multiply the equation by r and integrate over (0, R) to arrive at 

R{-u'{R)) + (N -2){u{0) -u{R)) ^ / rg{r)dr. 

Jo 

If now N = 2, then one has 

Ri-u'iR)) = / rgir)dr < "-^H^^^^ 



Dividing by R and integrating the result over (0, R) gives the claim. 
Now take iV > 3. Since —u'{R) > we can drop a term to arrive at 



(N - 2)(.(0) - uiR)) < / rgir)dr ^ ^ J^^ ^d. < ^ ^J^^ ^-^d. 



and then use Lemma 13.11 to evaluate the integral on the right and finish the proof. □ 

We now come to the result which will yield our upper estimates on u* . 

Theorem 3.3. Suppose u is a smooth semi-stable solution of {P\) on the unit ball B in M.^ , where 1 < 
N < 11. Then, for max{0, < t < 2 + ^/6, we have 



I 



' R^'-'dR < 1 (^ 4(2t + l) \^ 

Remark 3.2. Note that the above theorem only shows that is bounded away from 1 if 4i + 6 — 2N > 

2N — 1 which, once coupled with the other condition on t cannot be satisfied in the higher dimensions. This 
is to be expected since the extremal solution u* satisfies u*(0) = 1 for N > 8. 

Proof. Suppose u is a smooth semi-stable (so radial) solution of (Pa)- Then, the above linear estimate 
applied with g{r) := A(l — u)^^, gives that for all R E (0, 1), 

1 - UiR) < 1 - .(0) + >^^ir.Nmi-u)-^h^R^-f^ 



Now use the upper bound A* < "^^l^^^ from Proposition 12.11 take r = i + |, and replace the L"^— norm on 



the right via the energy estimate from Lemma 13.31 to obtain 

.X 4Ai(B)7(t+ |,iV) / 4(2i+ 1) \^ it+e-2^ 

i-.(i^)<i-.(o)+ UWi 
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This yields the inequality 

R^-UR R^-'dR 



1 - \\u\\l^ + 



4Ai(i3)7(t+j,JV) / 4(2t+l) ^ - pli±| 
27 \^4t+2-t2 J ^ 



2 s2t+3 — ^'"^JV I n — 1/1' R^l "124+3 ■ 

4Ai(i3)7(t+ j,JV) / 4(2t+l) \ " j^ ^nfZ?'^ \ "^^ ^ ^ -^-^ 



But the right hand side is actually equal to — u) || , hence we can again use the energy estimate 
from Lemma 13.31 to majorize it and complete the proof. 

□ 

Remark 3.3. Using Maple to approximate the integral in il9\) while optimizing over t, we get the following 
estimates on the extremal solution u* of {P\) on the unit ball in R^. 

1. If N = 1, then ||w*||l~ < -49... 

2. If N = 2, then ||w*||l- < -55... 

We now obtain some explicit upper bounds on u* in dimensions N = 1,2. For that, we define 

4Ai(S)7(t+|,iV) / 4(2i+l) 



C(t,N):= 

^ ' ' 27 V 4t + 2 - < ^ 

Corollary 3.1. Suppose u* is the extremal solution of {P\) on the unit ball in R^. 

1. IfN^ 1, then 

!f 2t + 3 \ 2j- + 2 

2. IfN = 2, then 

Proof. 1) For < i < 2 + \/6 one has > 1 and so we can replace the power on R in (fTOj) by 1, so as 

to be able to explicitly calculate the integral in p^ . One then drops a few positive terms to arrive at the 
desired result. 

2) For i < < < 2 + \/6 one has '^*2f+3^^ > 1, so again we replace the power on R in by 1 and carry 
on as in the first part. □ 



4 Effect of power-law profiles on pull-in distances 

Our goal in this section is to study the effect of power-like permittivity profiles f{x) = \x\" on the problem 
{P\.a) (our notation for (Pa,|x|<»)) on the unit ball B — -Bi(O). Numerical results - in particular those 
obtained by Guo, Pan and Ward in jTSl for MEMS nonlinearities- give lots of information, but the most 
intriguing one is their observation that on a 2-dimensional disc, the pull-in distance does not depend on a, 
at least in the case where F{u) = (1 — and that the solution develops a boundary-layer structure near 

the boundary of the domain as a is increased. In other words, the norm of the extremal solution of 
{M\^a) is independent of a > 0. In this section, we shall give a simple proof of this observation and other 
interesting phenomena, which actually holds true for more general nonlinearities. 

We first observe that since r — > r" is increasing, the moving plane method of Gidas, Ni and Nirenberg 
]T2\ does not guarantee the radial symmetry of all solutions to {S\j). However, one can show as in [TU] the 
following proposition. 
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Proposition 4.1. Let ^ be a radially symmetric domain and assume f is a radial profile on Q. Then, the 
minimal solutions of {P\j) on £7 are necessarily radially symmetric and consequently 

A*(ri, /) = A*(ri, /) = sup {A; {P\j) has a radial solution}. 

Moreover, if fl is a ball, then any radial solution of {Pxj) attains its maximum at 0. 

Proof. It is clear that A*(ri, /) < A*(ri, /), and the reverse will be proved if we establish that every minimal 
solution of (-Pa,/) with < A < A*(ri, /) is radially symmetric. The recursive linear scheme that is used to 
construct the minimal solutions, gives a radial function at each step, and the resulting limiting function is 
therefore radially symmetric. 

For a solution u{r) on the ball of radius R, we have Ur{0) = and 

iV - 1 

—Urr Ur — Xf{r)F{u) in (0, R) . 



u(r) attains its maximum at 0, and that - just as in the case where / = 1 - we have ||w*||oo — u*{0). □ 



Hence, -^^^-j^r^ = Xf{r)r'^-'^F{u) > 0, and therefore < in (0,i?) since Ur{0) = 0. This shows that 



It follows from this proposition that for radially symmetric domains O and profiles /, the extremal 
solution u* is necessarily radially symmetric and that the pull-in distance is nothing but u*(0). We shall 
denote by A* (TV) (resp., m*) the pull-in voltage (resp., the extremal solution) of (-Pa,/) when f{x) = 
and O is the unit ball in M.^ . 

We now make the following crucial observation. 

Proposition 4.2. For any a > ~2, the change of variable u{r) — ii;(r^+^) gives a correspondence between 
the radially symmetric solutions of the equation 

-Aatu^ A(l-Hf)2|x|"P(u) inB, 

u = on dB, ^ ' 

in dimension N and those of the equation 

I —/S.2(N+a,)W — \F{w) inB, ^ ^ 

w — Q on oB, 



the potentially fractional - dimension N(a) — ■ Moreover, we have 

A*(7V) = (l + f)2Aa(7V(a)) and uUr) = w* {r^+ ^ ) , (22) 
where m* is the extremal solution for i20\) and vu* is the extremal solution of \21\l . 

Proof: Indeed, by noting that for a radially symmetric u, we have Ajvw — u" + ^^^ u' , a straightforward 
calculation gives that 

^Nu{r) + (1 + ^fXr-F{u{r)) = (1 + ^fr^ (^i,^^)w{r^+^) + \F{w{r^+ 

where N{a) — ^^^q^^- The rest follows from the uniqueness of the extremal solutions. 

The above transformation allows us to deduce many results for the case of a power-law profile, from 
corresponding ones associated to constant profiles. The fact that it preserves the -L°°-norm has consequences 
on the pull-in distance and on the role of the profile in the critical dimension. It does also give proofs 
for various intriguing phenomena displayed by the numerical results below, especially in the case of a two 
dimensional disc, where the transformation does not alter the dimension since then N{a) — 2. 

The following corollary summarizes these consequences. 
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Corollary 4.1. With the above notations, the following hold: 
1. For any dimension N > 1, we have for a >> 1, 



A;(A^)~(l + ffAS(2). (23) 



2. IfN = 2, then 



A;(2) = (l + f)2Aj;(2) and \\u*J ^ \\u*J for all a > -2 . (24) 

Proof 1) From the above proposition, we have X*^{N) = (1 + f )^A5(2^±fii), and A5(2^±|i) - A5(2) 
whenever a is large. 

2) follows from the fact that for = 2, we then have Na = 2 for each a which means that 

A:(2) = ^^^AS(2), (25) 

and the pull-in distance in dimension 2 on the ball is ||l=o = ||w*||ioc, where u'^{r) = w*{r^~^^). The 
pull-in distance is therefore independent of a. □ 

Corollary 4.2. The following estimates hold in a MEMS model with a power-law permittivity profile, i.e., 
if F{u) = (1 - u)-^ and f{x) = |a;|". 

1. For any dimension N > 1, we have for a >> 1, 

A;(A^)- 0.789(1 + |)2. (26) 

2. IfN^ 2, then 

A;(2) = 0.789(1 -f f)2 and \\u*Jl^ = 0.U5 for all a > -2. (27) 

3. If 1 < N < 7 or if N > 8 and a > aN -.^ ^^4^2^^ ' ^'^^^ extremal solution u* of (Mx,a) on the 
hall is classical and the pull-in distance H^oo < 1. 

4. If the dimension N > 8, and < a < '.— '^^4~+2Vi^ ' then the extremal solution is exactly 
w* (x) = 1 — |a;|~5~^ which means that 

KiN) ^ and IKIUoo^l. (28) 

Proof. 1) and 2) follow from the above proposition and the fact that Aq(2) = 0.789 and ||uQ||ioo = 0.445. 

3) The extremal solution u* of (il^A.a) is regular if and only if HuallL^ = ||ii'*||L°° < 1, where w* is 

the extremal solution for {M\) in dimension N{a). According to [TD], this happens if < 1 _|_ | ^ 
which means that a > on := M^l±:^, 

4+2\/6 

4) Note first that u*(x) = 1 — |a;|T~ is a i7g(_B)— weak solution of {Mx^ix]") for any a > A — 3N. The 
voltage is then Aq(A^) = ^ Since now ||w*||ioc = 1, then by Proposition [Ol it remains only to 



show that for all (j> e II^{B), 



But Hardy's inequality gives for > 3 that Jg |V0p > ^"^^ Jg for any 4> G Hq{B), which means 
that ^ holds whenever 2XaiN) < L^L^ or, equivalently if > 8 and < a < oat = ^^^j^^f^f^- □ 
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The above scaling has also the following direct consequences. 

Corollary 4.3. Suppose F is a regular nonlinearity and < 10 + 4a, then the extremal functional u* of 
{Px.a) on the ball is classical. 

Proof. Cabre and Cappella [21 showed that the extremal solution on the ball is always bounded for iV < 9. 
They were only interested in integer dimensions, but an inspection of their proof indicates that the same 
result holds for any fractional dimensions < 10. Combining this with our observation in Proposition 14.21 
completes the proof. To see that this is optimal one recalls that when F{v) = e" the extremal solution 
is unbounded in N = 10. Using this fact and the change of variables above yields the optimality of this 
result. □ 

Remark 4.1. One can also use this change of variables to study permittivity profiles with negative powers 
(i.e., f{x) — \x\" for > a > —2. For example suppose F{u) = e", then using the above change of variables, 
one can show that for a fixed N (3 < TV < 9), the extremal solution associated with 

-Am = |a;|"e" on B, 

is singular for a £ (—2, ^°~^ ], while it is a classical solution for a € ( ^'^^^ , 0). 



5 Asymptotic behavior of stable solutions near the pull-in voltage 

We now establish pointwise upper and lower estimates on the minimal solutions ux in terms of A, A*, the 
extremal solution u* and ^^^a |_,^__,^. • For simplicity we restrict our attention to F(u) — e" and F{u) — 
(1 — In addition we allow fractional dimensions for results on the unit ball since then, one can apply 

the results of the previous section to deal with power-law profiles (Pa,q)- We first recall that by using 
Proposition ll.il one can show the following: 

• If F{u) = e", then u*{x) — log(j^) is an extremal solution on the unit ball in at A* = 2N — 4, 
provided > 10. 

• If F{u) = (^i\yi , then u*{x) — u*{x) — 1 — \x\i is an extremal solution on the imit ball in at 
A* = provided N > 

Theorem 5.1. Let u* denote the extremal solution of {P\) on a smooth bounded domain f2 in R^. 
1. If F{u) = (1 —u)^^, then for < A < A*, we have 

ux{x) < ( ) ' u*{x) for a.e. x e Q. (30) 



^A , 

Moreover, if is the unit ball in R^ with N > MMv^ ^ 7.93...^ then /or < A < A* = we have 

, , 2 3(A* — A) /, , -Will V9iV^-84W + 100 \ /A\^,^ i|2, ,„^, 

(6jv_8) -Ij <ux{x)<i^—j (l-|x|3), (31) 

for a.e. x € fl. 

2. IfF{u) = e", then for < \ < X* , 

uxix) < log j for a.e. X en. (32) 
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Moreover, if O is the unit ball in M with N > 10, then for < A < A* = 2A'' — 4 we have 



'^'^W^ - W^) [\^\-^'^^^ - ^ "a(-) ^ log , (33) 

for a.e. x Gfl. 

Proof. The upper estimates follow easily from the minimality of u\ and the fact that a; i— > (■^) u*(a;) (resp., 
X I— > log (^ )^,_^^^^-u* ^ is a supersolution of (Px) in the case that F{u) = (1 — (resp., F{u) = e"). 

For the lower bound, we shall proceed as follows: First, recall that A ua is diffcrentiable and increasing 
on (0, A*), and so if one defines vx := j^ux, then vx is positive and solves the linear equation 

-Av = F{ux) + XF'{ux)v in f], 
v = on dn, 

where, F is given by either e" or (1 — m)~^. We shall need the following notion. 

Definition 5.1. An extremal solution u* associated with {Fx) is said to be super-stable provided there exists 
e > such that 

(A* + e) < I for all V e H^in). 

Note that if u* is a super-stable extremal solution then /ii(A*, u*) > 0. We shall see at the end of this 
section that the converse is however not true. We first estabhsh the following result. 

Lemma 5.1. Assume ft is a smooth bounded domain in . Then, 

1. For < A < A*, Vx is the unique HQ — weak solution of (Qx)- 

2. X>-^ vx is increasing on (0, A*), and therefore v*{x) := lim>i^A* 'Vx{x) is defined for a.e. x G SI. 

3. A I— > ua is convex on (0, A*), and therefore for < A < A* we have for a.e. x G ft, 

ux{x)>u*{x) + {X-X*)v*{x). (34) 

4- If u* is super-stable, then v* is the unique H^^—weak solution of{Q)x'. 

Proof. (1) One can use the fact that /Ui(A,ma) > 0, and a standard minimization argument to show the 
existence of an ffg "Solution to {Qx}- Using the fact that /ii(A, ua) > one can see that the solution is 
unique. 

(2) Let < A < A* and e > small. Note first that 

-A{vx+e-Vx) = F{ux+e)-F{ux)+eF\ux+e)vx+e 
+ XF'{ux+e)vx+e " XF' {ux)vx 

= g{x) + XF'{ux){vx+e -Vx), 

where 

g{x) := F{ux+e) - F{ux) + sF'{ux+e)vx+e + X{F'{ux+e)vx+e - F'{ux)vx+e) 

is in H^{Q) and is positive. Now set w := vx+e — vx in such a way that w solves 

—Aw = g{x) + XF' {ux)w on O, 
w = on dO,. 

Testing this equation on w~ gives 

- / gw~ >iii{X,ux) / {w~f, 
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and hence w~ — a.e. in fi. By the maximum principle one then get that w > in f2 and hence that 
A — > is increasing. We can therefore define the hmit v*{x) := hm^^A* 'Vxix), which exists a.e. x in Q, 
though it might be infinite on a large set. 

(3) The convexity of A i-^ ua follows from the fact that A i-^ wa is increasing. We can therefore write 
u\ > Ut + {X — t)vt for < A, t < A* and a.e. x Q fl. The claim now follows by letting t go to A*. 

(4) Since u* is super-stable one has 

(A + £) / i^'(?iA)^' < / |VVP VV^ e Hi 
Jn Jn 

Using this and testing {Q\) on v\ gives 

e [ F'{ux)vl < [ F{ux)vx. 

Since F is either F{u) — e" or F{u) = (1 — m)^^, the left hand side is necessarily bounded. From this and 
again by testing (Qa) on v\ one sees that v\ is bounded in Hq. Passing to limits, one sees that v* is a 
Hq— weak solution of (Qa*)- The uniqueness follows from the fact that fii{X*,u*) > 0. □ 

We now complete the proof of Theorem 15.11 For that we assume that ft is the unit ball in M^. It is 
then easy to show using Hardy's inequality that the explicit extremal solutions for (Pa) given above, are 
super-stable provided TV > 10 (resp., N > ^ 7.93...) when F{u) = e" (resp., F{u) = (1 - uy^). An 

easy calculation also shows that 

(when F{u) = e") resp., 

v*(x)^ b ~+^+ 6 -1 , 

( when F{u) = (1 — w)^^) are iJp— weak solutions of (Q)a* in the respective cases, assuming the dimension 
restrictions above. Using this and the earlier convexity result gives the desired lower bounds for > 10 {N > 
14+4 -| ^Yie exponential and MEMS cases respectively. To obtain the result for the critical dimensions 
one passes to the limit in A^. We omit the details. □ 
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